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Abstract. We report the discovery of a new gravitationally 
lensed QSO, at a redshift z = 1.689, with four QSO com- 
ponents in a cross-shaped arrangement around a bright galaxy. 
The maximum separation between images is 2^6, enabhng a 
reHable decomposition of the system. Three of the QSO com- 
ponents have g 19.6, while component A is about 0.6 mag 
brighter. The four components have nearly identical colours, 
suggesting little if any dust extinction in the foreground galaxy. 
The lensing galaxy is prominent in the i band, weaker in r 
and not detected in g. Its spatial profile is that of an elliptical 
- - galaxy with a scale length of ~ 12kpc. Combining the mea- 

■ sured colours and a mass model for the lens, we estimate a 
' most likely redshift range of 0.3 < 2: < 0.4. Predicted time de- 

fS| lays between the components are ^ 10 days. The QSO shows 
evidence for variability, with total g band magnitudes of 17.89 
^ ^ and 17.71 for two epochs separated by ~ 2 months. However, 
_ the relative fluxes of the components did not change, indicating 

■ that the variations are intrinsic to the QSO rather than induced 
, by microlensing. 



fSJ ] Key words: Quasars: individual: HE 0435—1223 - Quasars: 
general - Gravitational lensing 
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1. Introduction 

The majority of known gravitationally lensed QSOs displays 
image splitting into two components. Such systems offer rel- 
atively few constraints for their mass distributions - usually 
just the two positions since the flux ratios might be changed 
by microlensing. Both for the purp ose of studying len sing po- 
tentials (Keeton, Kochanek & Falco Keeton et al.l998) and for 



the purpose of measuring cosmological parameters (Schechter 
Schechter2001), quadruply imaged quasars are considerably 
more useful than their doubly imaged counterparts. The dis- 
covery of a new quadruply imaged is therefore welcomed by 
everyone except those who struggle to explain why their rela- 
tive num bers are so high in radi o lensing surveys (Rusin and 
Tegmark [Rusin & Tegmark2()()l ) . 

In this paper we report the discovery of a new gravitation- 
ally lensed QSO with quadruple image splitting. The object 
was found to be multiple as part of a new, ongoing imaging sur- 
vey for lensed quasars using the Magellan Consortium's 6.5 m 
Baade telescope on Cerro Las Campanas. We present the first 
spectrum of this QSO and analyse gri imaging data to estab- 
lish a first suite of astrometric and photometric properties of 
the QSO components and the lensing galaxy. We then discuss 
the available constraints on the surface mass distribution in the 
lens and conclude with some prospects for future observations. 



2. Observations 

2.7. Identification and spectroscopy 

HE 0435 — 1223 was originally found as a high-probability 
QSO candidate, selected in the course of the Ham- 
burg/ESO digital o bjective prism survey (Wisotzki et al. 
Wisotzki et al.2000 ). Subsequent low-resolution spectroscopy 
was performed in March 2000 using the ESO 1.52 m telescope 
on La Silla with its Boiler & Chivens spectrograph. A 5 min 
exposure was enough to confirm the nature of this object as a 
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* Based in part on observations obtained with the Baade 6.5-m tele- 
scope of the Magellan Consortium. Also based in part on observations 
collected at the European Southern Observatory, La Silla, Chile. 



Fig. 1. Spatially unresolved slit spectrum of HE 0435— 1223, 
taken prior to the Magellan images in order to confirm the QSO 
nature of the source. 



bona fide QSO, and to determine its redshift z = 1.689 (mea- 
sured from the peak of the C IV emission line). The spectrum 
is displayed in Fig. [l] and reveals a fairly average QSO with- 
out obvious peculiarities. At that time, only a photographic sky 
survey image existed, showing the QSO as an unresolved sin- 
gle source. It is now clear that only a fraction of the quadruple 
image ensemble was located within the 2" spectrograph slit, 
and the spectrum therefore constitutes an ill-defined blend of 
all four components. However, the very similar optical colours 
(see below) make it likely that Fig. |l| is representative for each 
of the QSO components. 

2.2. Imaging 

High resolution optical images were obtained at the 6.5 m 
Baade telescope on 14 Dec 2001, equipped with the Magellan 
Instant Camera (MagIC) and a 2k x 2k CCD. The image scale 
at the f/1 1 Nasmyth focus was approximately 0'.'0692 per pixel. 
A Shack-Hartmann wavefront sensor (Schechter et al. 2002, in 
preparation) was used to update the focus and translation of 
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Fig. 2. Gallery of images of the quadruple QSO, with g, r and i bands, from top to bottom. The left two panels show the 
unprocessed image with different linear cut levels. The two right-hand panel show the deconvolved images, with and without the 
quasar components. Each image measures 8^8 x 8^8 in the sky; the arrows in the top left panel are 2'!^ long. Notice the possible 
companion to the lensing galaxy to the SW. 



the secondary and the twelve most flexible elastic modes of the 
primary mirror at half minute intervals. 

After a first r' band 'snapshot' of 30 sec which revealed the 
multiple nature of the QSO, five further images were obtained: 
One each in the SDSS g' and i' bands, and three in r^ of 300 sec 
exposure each (for simplicity, we shall refer to these bands as 
gri in the following). The effective seeing was between 0^6 and 
0^8 (for the g image). At the scale of 0V07/pixel, the images 
were highly oversampled. In Fig. ^ we show postage stamps 
of the QSO in the three bands. The symmetric configuration 
of four bright point sources plus a fuzzy object in the centre 
immediately recalls the famous 'Einstein Cross' Q 2237+0305 



(Huchra et al. Huchra et al.l985), and the image alone leaves 



little doubt that HE 0435 — 1223 is another case of a quadruple 
gravitational mirage. Given the nearly textbook arrangement 
of the components in HE 0435—1223, we shall assume hence- 
forth that there is no reasonable alternative to the lensing hy- 
pothesis in this object. A colour composite image is shown in 
Fig. I 



Table 1. Differential astrometry for HE 0435-1223. Compo- 
nents A-D are labelled clockwise, starting with the brightest 
component A. 



Component 








[arcsec] 


[arcsec] 


A 


0.000 


0.000 


B 


-1.483 ± 0.002 


0.567 ± 0.002 


C 


-2.488 ± 0.003 


-0.589 ± 0.002 


D 


-0.951 ± 0.001 


-1.620 ± 0.001 


G 


-1.15 ± 0.05 


-0.51 ± 0.05 



Second epoch data was secured on 15 Feb 2002, using the 
same telescope and instrument configuration as above. Three 
120 sec images in the g band were taken at 1" seeing under pho- 
tometric conditions. The standard star sequence PG 1047+003 
was observed to provide a zeropoint for the g band. 
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i' and r' bands. The point sources were over sub traded, and 
the source positions were systematically shifted inwards, both 
clearly caused by the fact that the contribution of the lensing 
galaxy has not been taken into account in the fitting. 

We have therefore chosen to apply a more sophisticated 
method. The algori thm devised by Magain, Courbin & Sohy 
( Magain et al.l998 i hereafter MCS) combines multiple point 
source fitting with a locally regularised deconvolution of the 
smooth image component. The algorithm has been demon- 
strated to yield astrometrically and photometrically u nbiased 
results (MCS; see also Burud et al. [Burud et al.l998 ) and is 
thus well- suited to our case. We have used an 'alpha release 
public version' of the code to decompose each image into four 
point sources plus a smooth 'background' channel. The recon- 
structed images are shown in the right-hand panels of Fig. ^ 
adopting a final resolution of 0'.'21 FWHM in r and i, and of 
0'.'28 in The lensing galaxy is very prominent in the recon- 
structed i band image, fainter but clearly detectable in r, but 
swamped by PSF residuals in g. The non-detection in g is con- 
firmed by the DAOPHOT analysis which also shows no signif- 
icant galaxy residual at the expected location. 

While the MCS method is probably close to optimal in its 
treatment of the point sources, we found that the exact shape of 
the reconstructed lensing galaxy depended critically on the lo- 
cal smoothing term A which basically is a free parameter in the 
algorithm; slight variations of A (by a factor of ~ 2 or so) led 
to significantly different scale lengths. We therefore reverted to 
a hybrid approach: Using the results of the MCS routine for the 
point source positions and fluxes, we went back to DAOPHOT 
and performed another PSF subtraction cycle. The resulting im- 
age of the lensing galaxy looks almost free of PSF residuals and 
certainly good enough to perform a first reliable morphological 
analysis of the lensing galaxy. 



Fig. 3. Colour composite images of the unprocessed (left) and 
MCS-deconvolved (right) Magellan data. The RGB colour 
channels are represented by the i, r, and g band data, respec- 
tively. Orientation as in Fig. |[ 



3. Analysis 

3.1. PSF fitting and deconvolution 

In order to decompose the observed images into four point 
sources plus an extended contribution from the lensing galaxy, 
we followed two independent approaches. First, we performed 
a straightforward PSF fitting analysis using DAOPHOT II 



(Stetson et al. Stetsonl987). This confirmed the visual impres- 



sion that components A-D (numbered clockwise starting from 
the brightest component) are indeed point-like, but the PSF- 
subtracted image looked less than satisfactory especially in the 



3.2. Astrometry 

Rectangular positions of the four point source components 
were provided by the MCS PSF fitting/deconvolution code. Ar- 
bitrarily adopting the brightest component A as reference point, 
we Hst the relative positions of HE 0435-1 223 B-D in Table |l|. 
Errors were estimated from the rms deviations between individ- 
ual frames as well as between individual deconvolution runs, 
which differ by adopting slightly different parameters for the 
smoothing scale. All measurements agree with each other ex- 
traordinarily well, with deviations from the mean of never more 
than a few milli-arcseconds, with no systematic effects visible 
between the different photometric bands. 

The lensing galaxy had to be treated separately, as it is not 
modelled as an individual 'object' by the MCS code. However, 
in the r' and i' bands it has a well-defined centroid which we 
measured by fitting a Gaussian to the inner region. When ex- 
pressed relative to image A, the galaxy centroid is shifted by 
~ 0V05 between i and r. We estimate the uncertainty of the 
galaxy position to be at least of this order. 

In order to relate the above relative astrometry to a com- 
monly used reference system, we have tied the measurements 
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Table 2. Differential photometry for the four point source com- 
ponents, expressed as magnitudes relative to the total aperture 
flux in each band. The g band values are equally valid for first 
and second epoch data. 



Component 


9 


r 


i 


A 


1.10 


1.10 


1.11 


B 


1.61 


1.61 


1.59 


C 


1.70 


1.67 


1.63 


D 


1.74 


1.78 


1.82 



Table 3. Photometry for HE 0435-1223 (first epoch). CaH- 
brated magnitudes of the QSO (including galaxy) and its com- 
ponents are relative to the reference star as described in the 
text. 



Component 


9 


r 


i 


9-r 


r — i 


QSO total 


17.89 


17.15 


16.53 


0.74 


0.62 


A 


19.00 


18.44 


17.95 


0.56 


0.49 


B 


19.51 


18.95 


18.43 


0.56 


0.52 


C 


19.60 


19.01 


18.47 


0.59 


0.53 


D 


19.64 


19.10 


18.66 


0.54 


0.44 


G(oo) 




19.16 


18.05 




1.11 


G(0'.'7) 


>22.70 


20.91 


19.82 


>1.80 


1.08 


Ref. star 


17.80 


16.56 


15.98 


1.24 


0.58 



to the Digitized Sky Survey (DSS) astrometric solution. The 
quasar itself is unresolved in the DSS images, but four nearby 
stars are found in all our CCD images as well as in the DSS. 
We computed simple bilinear solutions for each CCD image 
and obtained the following coordinates of Component A: R.A. 
= 04^ 38"^ 14!90, Dec = -12° 17' 14'.'4 (J2000.0), with errors 
defined entirely by the unknown systematics of the DSS solu- 
tion in this area. 



3.3. Photometry 

Relative photometry of components A-D is directly provided 
on output by the MCS code. Even without any calibration it 
is clear that all components have very similar colours. In Ta- 
ble H we present the fractional contributions of each compo- 
nent to the total QSO flux (without lensing galaxy), expressed 
in magnitudes. Between the photometric bands, the contribu- 
tions of components A and B show deviations of no more 
than ±0.01 mag, while components C and D show marginally 
significant deviations but are still very nearly independent of 
wavelength, within ±0.05 mag. Note that the numbers for the 
g band are identical between first and second epoch to within 
±0.01 mag. 

For measuring the brightness of the lensing galaxy we con- 
structed PSF- subtracted images with DAOPHOT, but using the 
relative astrometry of the QSO components from the MCS de- 
convolution. The resulting galaxy-only images look quite clean 
and devoid of strong PSF subtraction residuals, except for the g 



band where we regard the lensing galaxy to be undetected, be- 
cause of the significant PSF residuals in that image. We quote 
two values for the galaxy flux: one obtained by integrating to 
infinity and a second through a small aperture of 10 pixels ra- 
dius (0''7). The small aperture was included because within this 
radius, residual PSF contamination should be close to negligi- 
ble. While a significant fraction of the total galaxy flux will be 
missed this way, the galaxy colours should be rather reliable. 
For the ^ — r colour, the formally measured g band aperture 
flux of the galaxy is treated as an upper limit. 

Since the first epoch data were obtained under nonphoto- 
metric conditions, we have tied the g band photometry of the 
QSO to a nearby star, visible within the CCD field of view and 
located at R.A. = 04^ 38"^ 12^9, Dec = -12° 17' 52'' (J2000.0). 
We found this star to have g = 17.80 in the second epoch 
data. Assuming the star to be non- variable, the total (large aper- 
ture) magnitude of the QSO is ^ = 17.89 a the first epoch and 
g = 17.71 at the second epoch, with uncertainties of at most 
±0.02 mag. There is thus evidence for significant variability of 
the QSO. However, as stated above, the fractional contributions 
of the individual QSO components to the total g band flux are 
identical in first and second epoch data. We return to this issue 
below. 

As the second epoch provided only calibrated g, but no 
r and i data, we have worked out at least approximate ze- 
ropoints for these latter bands as follows. We found that the 
above reference star has an entry in the database of digital ob- 
jec tive prism spectra of the Hamburg/ES O Survey (Wisotzki e t 
al. [Wisotzki et al.2Q0Q| ). ChristHeb et al. ( [ChristHeb et al.200l| ) 
showed that these spectra can be used to obtain (B — V) colours 
of normal stars with an accuracy of better than 0. 1 mag. Follow- 
ing their recipe, we estimate 5 — F 1 . 1 for our star and then 
used a matching model atmosphere spectrum[] to estimate its 
colours in the SDSS system to ^ — r ^ 1.24, and g — i c:^ 1.82. 
Combining the second epoch g band calibration with the es- 
timated intrinsic colours of the reference star, we could thus 
tie the first epoch r and i band data to an external photometric 
system. The uncertainty of this 'calibration' procedure is con- 
siderable, leading to zeropoint errors of ~ 0.15 mag for r and 
~ 0.2 mag for i (assuming a standard error of ±0.1 mag in 
B — V of the reference star). However, the errors in r and i are 
strongly correlated, and the zeropoint uncertainty in r — i colour 
is only ~ 0.06 mag. Results for the final calibrated photometric 
results are collated in Table |[ 



4. Discussion 

4.1. Properties of the lensing galaxy 

Lacking direct spectroscopic evidence, we have estimated the 
redshift of the lensing galaxy using its measured colours. Fig- 
ure Q shows the redshift-colour dependence of three simple 
galaxy models, represented as single stellar populations (SSPs) 



^ obtained from R. Kurucz' web site at 

http : / /cf aku5 . harvard . edu/ 
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Fig. 4. r — i vs. g — r colour variations of single stellar popu- 
lation models, as a function of redshift. Steps of Az = 0.1 are 
denoted by the small tickmarks (filled symbols: population age 
14 Gyrs; open circles: 4 Gyrs; open squares: 1 Gyr). The ar- 
row+bar symbol represents the measured r — i and upper limit 
on ^ — r of the lensing galaxy. The inclination of the error bar 
indicates the correlated errors in the two colours. 




1 2 

radius [arcsec] 



Fig. 5. Radial surface brightness profile of the lensing galaxy, 
measured in the PSF- subtracted i band image. The solid line 
shows the best-fit de Vaucouleurs profile (inner regions extrap- 
olated as shown by the dotted line). 



side the small (0'!7) aperture, with the width and inclination of 
the bar reflecting chiefly the uncertain calibration (recall that 
the errors inr — i and g — r are correlated). 

This plot suggests that the redshift is probably around z = 
0.4-0.5, the best estimate depending on the unknown intrinsic 
spectrum of the galaxy. The youngest (1 Gyr) SSP model yields 
the highest redshift (zi = 0.55) and is only marginally consis- 
tent with our upper limit for the ^-band flux. With the more 
reasonable assumption of an old stellar population dominating 
the SED of the lensing galaxy, the best redshift estimate drops 
to z = 0.41 for a 14 Gyrs old SSP. Because the 14 Gyrs track 
bends sharply in exactly this part of the colour-colour diagram, 
the corresponding formal error is highly asymmetric. Adopt- 
ing a 2(7 confidence limit, we predict a redshift in the interval 
0.3 < zi < 0.5; other arguments detailed below will show 
that it is more likely to be near the lower end of that range. Of 
course, two colours are insufficient to distinguish an evolved 
population from a younger one that involves dust extinction, 
but a very dusty lensing galaxy can probably be ruled out from 
the fact that the four QSO images have such similar colours. 

For the formally best redshift estimate of zi = 0.41, the 
absolute magnitude My of the lensing galaxy can be converted 
directly from the i band measurement without any K correc- 
tion. We obtain My = -23.8 for Hq = 65 km s"^ Mpc-\ 
Qrn = 0.3, ft A = 0.7, about a factor of 5 brighter than an 
field galaxy. This would drop to Mr = —23.0 or ~ 2L^ for 
zi = 0.3. 

The reconstructed i band image, without the point source 
components, allows us to investigate the gross morphological 
properties of the lens. Figure || shows the radial surface bright- 
ness profile, measured in concentric annulli of 2 pixels width 
(0V138). Inside of ~ 0%, the profile is blurred by the seeing, 
but outside this radius it is clear that the galaxy follows approx- 
imately a de Vaucouleurs r^/^ law. The best-fit scale length is 
Te = 2'.'0 corresponding to (12 ± 3) kpc at 2;^ = 0.41 (red- 
shift uncertainty not included in error estimate). Notice that 
the galaxy is nearly perfectly round, with isophotes having 
1 — 6/a<0.1 everywhere. 

A possible companion is found 5'' to the SW of the lensing 
galaxy, at a position angle of 217°. It is visible in all bands, 
although the detection in g is only marginal. Its colours, mea- 
sured over an aperture of radius 10 pixels or 0'.'7, are r — z :^ 
1.4 ± 0.1 and g — r 1.6±0.2 (again with strongly correlated 
errors). The colours are thus consistent with the 'companion' 
being at the same redshift as the galaxy, but without spectro- 
scopic confirmation this conclusion must remain speculative. 
The total magnitude of the companion is i = 19.6, correspond- 
ing to Mv = -22.3 or - for z = 0.41. 



with ages of 1, 4, and 14 Gyrs, respectively. (For the lim- 
ited redshift baseline covered, we ignored passive evolution.) 
We used the recent sol ar metallicity models from Jimenez et 
al. ( limenez et al.2QQC| ), kindly provided by R. Jimenez. The 
bar+arrow symbol denotes the measurement of the galaxy in- 



4.2. The lensing potential 
4.2.1. Qualitative features 

The image configuration is roughly diamond shaped (indicat- 
ing perhaps a diamond in the rough), and nearly bilaterally 
symmetric with respect to the AC diameter. The AC diame- 
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Table 4. Observed and predicted quantities for the images and 
the lens. Positions are in arcseconds, relative to the observed 
lens centre, /i is the magnification factor, ti is a quantity related 
to the time delay (see text. Sect. 4.2.3). 



observed 
Aa A6 



predicted 
Aa A^ 



A 


1.150 


0.510 


1.150 


0.509 


7.4 


-0.798 


B 


-0.333 


1.077 


-0.334 


1.077 


-8.2 


-0.725 


C 


-1.338 


-0.079 


-1.339 


-0.078 


7.7 


-0.790 


D 


0.199 


-1.110 


0.201 


-1.110 


-5.0 


-0.636 


G 


0.000 


0.000 


-0.036 


-0.047 







ter is a factor of 1 . 14 times larger than the BD diameter. Taking 
the potential of the lensing galaxy to be roughly isothermal, 
this implies that the source of the quadrupole moment, whether 
due to a bar or a tidal pertur ber, lies along t he BD diameter 
(e.g. Witt, Mao and Schechter |Witt et al. 19951 ). The strength of 
the shear should be roughly half (AC - BD)/(AC + BD), or 
~ 0.07. The strength of the isothermal sphere model ought to 
be half the average of the two diameters, or 1.203 arcseconds. 



4.2.2. Quantitative aspects 

A model which works quite well for many lenses is a singular 
isothermal sphere with an external shear, as might be generated 
by the tide from a neighboring galaxy or clu ster. The two di- 
mensional projected potential (e.g. Kochanek Kochanekl991 ) 
is given by 

V^2d(6>) = 6r + |r2cos((/)-(/)^) (1) 

where b is the diameter of the isothermal sphere in arcseconds, 
r and (j) are the radial and angular parts, respectively, of the 
angular position 6 on the sky with respect to the galaxy, and 
(pj is the orientation of the shear on the sky, measured E of N. 
Given our sign convention, a bar or tidal perturber would be at 
right angles to (j)^ . In addition to the three explicit parameters, 
the position of the source adds two more parameters. 

Fitting this 5 parameter model to the positions (but not the 
fluxes) of the four images we find that the model gives rms 
errors of ~ 0'.'05. The rms residuals improve considerably if 
we take the position of the lensing galaxy to be free, adding 
two additional parameters. The predicted galaxy position is 
then some ~ 0'.'06 from the observed position, consistent with 
our measurement errors. We find b — 1'.'209, 7 = 0.074 and 

= 76.5°, with a source position (Aa, A(5) = (-0'.'0423, 
— 0'.'0076). The source of the tide is predicted to lie at P.A. 
-13.5° or 166.5°. 

In Table Q we give the observed and predicted positions 
for the four images, as well as the predicted magnifications, /i. 
Negative magnifications indicate the parity flip associated with 
saddlepoints of the time delay function. Components A and D 
are both observed to be brighter than predicted by ~ 0.5 mag 
relative to components B and C. 



For an isothermal sphere, the lens strength (measured in 
radians) is given by 

, Bis 4^^' 



(2) 



(Narayan and Bartelmann [Narayan & Bartelmannl999| ) where 
(J is the line of sight velocity dispersion of the isothermal 
sphere and where D^s and Dis are angular diameter distances, 
from observer to source and lens to source, respectively. Taking 
zi = 0.41 we get cr = 252 km/s. The Faber- Jackson relation 
would predict a substantially larger dispersion for a 5L^ galaxy. 
K\ zi = 0.30, we get a = 237 km/s, more nearly in line with 
the expectation for a 2L^ galaxy. 

The shear, 7 = 0.074, while not especially large for a 
quadruple lens, would call for an E3-E4 mass distribution if the 
flattening of the gal axy were the sourc e of the quadrupole mo- 
ment (Keeton et al. Keeton et al.l998 ). As the flattening of the 
light is unmeasurably small, this suggests that the quadrupole 
moment is indeed due to an external tide. In several quadru- 
ple systems the source of the tide, a cluster or group of galax- 
ies, is clearly vis ible at right a ngles to (Kundic et al. 
Kundic et al.l997|; Sch echter et al. [Schechter et al.l997| ; Kneib 
et al. Kneib et al.l998). In the present case the nearest galaxy, 
at P.A. 217°, is not at the expected position angle. There are 
several more galaxies in the neighborhood, again with no obvi- 
ous concentrations at the expected position angle. 

4.2.3. Time delay 

The gravitational time delay function, which gives the addi- 
tional travel time for an image at position Qi, is given by 

c Dis (206265)2 ^ ^ 

where (3 is the position of the source, and where we have ex- 
plicitly included the transformation from arcseconds to radi- 
ans. As the redshift of the lens is poorly known, in Table ^ we 
give only the square bracketted quantity, labelled ti, which has 
units of (arcsecond)^. We see that components A and C lead, 
followed by saddlepoints B and D. The longest delay, td — ta 
is 11?8 for zi = 0.41 or 7^9 for zi = 0.30. Predicted time 
delays are quite sensitive to the position of the lensing galaxy 
for symmetric lenses. The time delays predicted using the ob- 
served position for the lensing galaxy are 15 % smaller. 

4.3. Variability 



In Sect. 3.3 we demonstrated that the total QSO flux in the g 
band varied by almost 20 % within a timespan of 2 months. A 
comparison between our 2001/2002 CCD photometry and the 
original photographic discovery data reveals additional, albeit 
circumstantial, evidence for variability: On the UKST image of 
the Digitized Sky Survey, calibrated by photometric sequences 
taken in the course of the Hamburg/ESO Survey, we measure 
the QSO to have Bj = 16.2 ^ g Sit epoch 1984.9. On the 
spectral discovery plate obtained with the ESO Schmidt tele- 
scope (epoch 1996.9), the object has only Bj = 17.2, cali- 
brated against the same photometric sequence. While the HES 
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photographic photometry is well-defined only for point sources 
(with a global rms accuracy of ~ 0.15 mag), which might lead 
to a photometric bias for multiply imaged QSOs, inspection 
of the scan data revealed the quasar to appear perfectly unre- 
solved. We conclude that the differences between photographic 
and the more recent CCD measurements are most probably sig- 
nificant. 

For the two Magellan epochs Dec 2001 and Feb 2002, we 
can state that the variability is almost certainly intrinsic to the 
QSO and not caused by microlensing. This is given by the fact 
that the flux ratios of the four components have remained un- 
changed to within 1 %, despite an almost 20 % variation of to- 
tal flux. Microlensing-induced flux amplification would always 
affect each image differently, which is in contradiction to the 
observations. On the other hand, the short expected time delay 
means that intrinsic variations (assumed to be much slower) 
will show up quasi- simultaneously in all four images. 

While we can exclude microlensing as the cause of vari- 
ability between Dec 2001 and Feb 2002, this cannot be said 
about the origin of the much larger variations between 1986 
and today. In fact, the amplitude of more than 1 mag might 
even favour microlensing, as such amplitudes are rare for in- 
trinsic variation s of normal, n on-OVV QSOs (H ook et al. 
Hook et all 994 Cristiani et al. [Cristiani et al.l996| ). From the 
ex tensive monitorin g campaigns of the Einstein Cross ( Irwin et 
al. Ilrwin et al.l9"89i Wozniak et al. [Wozniak et al.2000t OGLE 
web site), on the other hand, it is known that microlensing am- 
plifications of more than one magnitude are possible and in fact 
happening. 

5. Conclusion 

The new quadruple QSO HE 0435 — 1223 is an almost textbook 
example for gravitational lensing, with its four nearly identical 
components arranged symmetrically around a luminous early- 
type galaxy. Unlike most other known quadruple systems, pho- 
tometric monitoring of this object should be relatively easy 
even in mediocre seeing conditions, because of its wide im- 
age separations. Furthermore, its location in the sky makes it 
accessible to both Northern and Southern observatories. 

Owing to its symmetry, the time delay is expected to be 
short, and accurate measurement of differential time delays 
might therefore be difficult unless the QSO should prove to 
be variable on very short timescales. This could limit the use- 
fulness of the object for cosmological purposes, but at the same 
time it makes it an attractive target for microlensing studies, be- 
cause of the relative ease to separate intrinsic and microlensing- 
induced variations. Notice that compared to the Einstein Cross 
Q 2237+0305, the higher lens redshift in HE 0435-1223 im- 
plies a ~ 10 X lower projected transverse velocity and hence 
a much longer characteristic timescale for high-amplification 
event from microlensing It will therefore be easier to obtain a 
well-sampled lightcurve, but unfortunately events will be rarer 
and take much longer to get covered. 

We have presented evidence that the QSO experiences sub- 
stantial flux variations on time scales of months and years. 



Whether microlensing could have a contribution in these vari- 
ations is not yet clear, but it can already be said with certainty 
that monitoring of HE 0435 — 1223 will be a promising task. 
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